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O ! ABSTRACT 

Q 

■ We report the discovery of two doubly-imaged quasars, SDSS J100128. 61+502756. 9 

and SDSS J120629.65+433217.6, at redshifts of 1.838 and 1.789 and with image sep- 
{Sj . arations of 2'.'86 and 2'.'90, respectively. The objects were selected as lens candidates 

from the Sloan Digital Sky Survey (SDSS). Based on the identical nature of the spectra 
I/*-) ' of the two quasars in each pair and the identification of the lens galaxies, we conclude 
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that the objects are gravitational lenses. The lenses are complicated; in both systems 
there are several galaxies in the fields very close to the quasars, in addition to the lens 
galaxies themselves. The lens modeling implies that these nearby galaxies contribute 
significantly to the lens potentials. On larger scales, we have detected an enhancement 
in the galaxy density near SDSS J100128.61+502756.9. The number of lenses with 
image separation of ~ 3" in the SDSS already exceeds the prediction of simple theoreti- 
cal models based on the standard Lambda-dominated cosmology and observed velocity 
function of galaxies. 

Subject headings: cosmology: observation — cosmology: theory — gravitational lens- 
ing — quasars: individual (SDSS J100128. 61+502756.9) — quasars: individual (SDSS 
J120629.65+433217.6) 

1. Introduction 

In strong gravitational lensing of quasars, the separations between multiple images, 9, is the 
most important observable linking observations to theory. Since the image separation is determined 
by the potential depth of the lens, the image separation distribution of lensed quasars offers a direct 
probe of the hierarchical structure of the universe. For instance, normal galaxies can produce 
strongly lensed quasars with image separations of ~ 1", while lenses with image separation > 10" 
can only be caused by clusters of galaxies. About 70 of ~ 1" lenses are known to date 1 , and there 
is one example of a lensed quasar system in which the lens potential is dominated by that of dark 
matter (Inada et al. 2003b; Oguri et al. 2004a). 

Among gravitationally lensed quasars, those with intermediate image separations (3" < 9 < 7" 
) are of great interest because they represent a link between small- and large-separation lenses. In 
the standard modeling procedure used to predict the distribution of image separations, assuming 
isothermal profiles and an a priori velocity function of galaxies, lenses with image separations > 3" 
are very rare, because even the largest early type galaxies do not have Einstein radii this large. 
Thus the probability for ~ 7" lensing is almost negligible. However, while Q0957+561 (9 = 6'.'26; 
Walsh, Carswell, & Weymann 1979) is primarily lensed by a galaxy, the image separation is boosted 
by the cluster in which the lensing galaxy resides. This example implies that the environment of 
the lens galaxy may significantly affect the distribution of image separations in the 3" — 7" range 
(Keeton, Christlein, &; Zabludoff 2000; Martel, Premadi, & Matzner 2002). In addition, a secondary 
mass along the line of sight could affect strong lensing (Wambsganss, Bode, &i Ostriker 2005), and 
this also may enhance the lensing probabilities in this image separation range. Finally, there is a 
predicted contribution in this range from clusters; simple theoretical models that include transition 



X A summary of known lensed quasar systems is available on the CASTLES homepage (Kochanek, C. S., Falco, 
E. E., Impey, C, Lehar, J., McLeod, B., & Rix, H.-W., http://cfa-www.harvard.edu/castles/) 
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of the property of lenses at ~ 1O 13 M (e.g., Oguri 2002) predict that large-separation lenses due to 
clusters begin to dominate the total lensing probability. Therefore, the overall lensing probability 
distribution for 8 > 3" is predicted to depend on the interplay of these two effects; the environmental 
effects and the emergence of cluster lenses. However, the overall lensing probability at 9 > 3" is quite 
small, thus a large number of quasars is needed to investigate the lensing probability distribution. 
Indeed, even the current largest homogeneous sample of lensed quasars (Myers et al. 2003; Browne 
et al. 2003) contains only one lens in this image separation range. 

In this paper, we present the discovery of two ~ 3" gravitationally lensed quasars, SDSS 
J100128.61+502756.9 (hereafter SDSS J1001+5027) and SDSS SDSS J120629.65+433217.6 (here- 
after SDSS J1206+4332). These gravitational lenses were identified from an ongoing lens search 
using the data of the Sloan Digital Sky Survey (SDSS; York et al. 2000; Stoughton et al. 2002; 
Abazajian et al. 2003, 2004, 2005). Currently the SDSS contains more than 50,000 spectroscopically 
classified quasars; thus the SDSS provides the opportunity to construct the largest homogeneous 
lens catalog in existence. Indeed, 13 new gravitationally lensed quasars have been found by using 
the SDSS (e.g., Inada et al. 2003a). In this paper, we describe photometric and spectroscopic 
observations of two new lens candidates and show that they are gravitational lenses. We model 
the lens systems and discuss the environments of the lens galaxies. We also compare the image 
separation distributions of lensed quasars in the SDSS (though still very preliminary because of the 
limited statistics) with a simple theoretical model. 

This paper is organized as follows. In §2, we briefly describe our method of searching for 
lens candidates from the SDSS data. Section 3 presents the results of both photometric and 
spectroscopic follow-up observations, and §4 shows the result of lens modeling. Section 5 is devoted 
to a discussion of the environments of the lens galaxies. We also discuss the lensing probability 
distribution, which is shown in §6. We summarize our results in §7. 

2. Selecting Candidates from the SDSS data 

All gravitational lenses presented in this paper were selected as lens candidates from the SDSS, 
which is a survey to image 10 4 deg 2 of the sky. The SDSS also conducts spectroscopy of galaxies and 
quasars that are selected from the imaging data (Eisenstein et al. 2001; Strauss et al. 2002; Richards 
et al. 2002; Blanton et al. 2003). A dedicated 2.5-meter telescope at Apache Point Observatory 
(APO) is equipped with a multi-CCD camera (Gunn et al. 1998) with five optical broad bands 
centered at 3551, 4686, 6166, 7480, and 8932 A (Fukugita et al. 1996; Stoughton et al. 2002). The 
imaging data are automatically reduced by the photometric pipeline (Lupton et al. 2001). The 
astrometric positions are accurate to about O'.'l for sources brighter than r = 20.5 (Pier et al. 
2003), and the photometric errors are typically less than 0.03 magnitude (Hogg et al. 2001; Smith 
et al. 2002; Ivezic et al. 2004). The spectra cover 3800-9200 A at a resolution of 1800-2100. 

We use spectroscopically classified quasars with z > 0.6 to search for gravitational lens candi- 
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dates. SDSS J1001+5027 and SDSS J1206+4332 are identified as lens candidates by our standard 
candidate selection algorithm (N. Inada et al., in preparation). This algorithm is based on the idea 
that the image of a quasar pair with a small separation appears to be more extended than that of 
single quasar, and characterizes the extent by the following SDSS image parameters: dev_L (the 
likelihood that the image of the object is fit by a de Vaucouleurs profile), exp_L (the likelihood by 
an exponential disk), and star_L (the likelihood by the point spread function). This algorithm has 
already found six new SDSS lenses (Inada et al. 2003a, 2005a,b; Pindor et al. 2004, 2005; Oguri et 
al. 2004b) as well as all previously known gravitational lenses in the SDSS footprint. However, the 
possible second lensed components of the candidates we study in this paper were also recognized 
as separate astronomical objects in the SDSS imaging data because of the relatively large image 
separations (Pindor et al. 2003). 

Figure 1 shows the SDSS i-band images of the candidates. In both images, the seeing was 
~ 1"2, better than the typical seeing of SDSS images, l'/35 (Abazajian et al. 2003). The Figure 
clearly shows two stellar components in each case; we designate the brighter ones as component A 
and the fainter ones as component B. Table 1 summarizes the result of the SDSS photometry of 
the components. In each system the components have colors similar to each other and the colors 
are consistent with that expected from quasars; this makes them excellent lens candidates. 

3. Follow-up Observations 

Follow-up spectroscopic observations were conducted on 2003 November 20 and 2004 June 21 
with the Double Imaging Spectrograph (DIS) on the Astrophysical Research Consortium (ARC) 
3.5-m telescope at APO. We used a 1"5 slit in low (blue) or medium (red) resolution mode. The 
resolution is about 2000. The final spectra cover the wavelength range of 3900 A to 9400 A. We 
aligned the slit along the direction joining the two images; since the separations between the 
images are large compared with the seeing size, extraction of the spectrum of each component was 
straightforward. The data were reduced using standard IRAF 2 reduction procedures. 

We also obtained deep images of the systems with the 8k mosaic CCD camera of the University 
of Hawaii 2.2-meter telescope (UH88) at Mauna Kea and with the Seaver Prototype Imaging Camera 
(SPIcam) of the ARC 3.5-m telescope. The UH88 imaging was conducted on 2004 May 23 and 25. 
We obtained VRI-b&nd images in the photometric nights; the exposure times were 360 seconds 
for each candidate and each band. The image scale was / /232 pixel -1 . We estimated the seeing 
as ~ Of! 8, better than the SDSS seeing size of ~ 1"2. The ARC z-band imaging was conducted 
on 2004 January 13. The exposure was 300 seconds, the seeing was ~ l'.'O, and the image scale 
was 0"282 pixel" 1 . Each frame was bias-subtracted and flat-field corrected. The magnitudes in the 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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UH88 images were calibrated by the standard star PG 1528+062 (Landolt 1992). Astrometry and 
photometry of the UH88 images are summarized in Table 2. 

3.1. SDSS J1001+5027 

Spectra and images of this object are shown in Figures 2 and 3, respectively. In Figure 2, 
both components show C IV, C III], and Mg II emission lines redshifted by z = 1.838. It is worth 
noting that the flux ratios (B/A) of the emission lines, particularly the C IV emission lines, are 
larger than that of the continuum; the equivalent widths of C IV emission lines are 56 A and 76 A 
for A and B, respectively. Such difference can be caused by the difference of the emission regions 
of continuum and broad emission lines combined with microlensing by stars. We find two extended 
objects near the quasar images (see Figure 3; the image separation is 2 / .'86). When subtracting a 
quasar component, we used a nearby star as a point-spread function (PSF) template. The objects, 
denoted by Gl and G2, have colors consistent with those of early-type galaxies at 0.2 < z < 0.5 
(Fukugita, Shimasaku, & Ichikawa 1995). Since galaxy Gl is nearly colinear with the two quasar 
components, it is likely that Gl is the main contributor to the lens potential. However, galaxy G2 
is also quite close to the lens system, and could affect the lens potential significantly. We note that 
component B could be reddened by the lens galaxy Gl. 

3.2. SDSS J1206+4332 

The spectra of the two components shown in Figure 4 show Si IV, C IV, C III], and Mg II 
emission lines at the same wavelengths, which supports the idea that this is a gravitationally lensed 
quasar at z = 1.789. However, the flux ratios (B/A) of the emission lines are slightly smaller 
than that of the continuum; the equivalent widths of C IV emission lines are 69 A and 53 A for A 
and B, respectively. The images in Figure 5 clearly reveal the lensing galaxy Gl as well as two 
quasar components separated by 2'/90, further supporting the lensing hypothesis. We also found 
other galaxies G2 and G3 near component B. These galaxies, particularly G2, may contribute to 
the lens potential to some extent. Indeed, both Gl and G2 are quite red (see Table 2) and are 
consistent with being high-redshift (z > 0.7) early-type galaxies (Fukugita et al. 1995), while G3 
is blue and thus may be a chance superposition of a local galaxy. There is strong Mg II absorption 
(equivalent width > 2 A) at ~ 4900 A in the spectrum of component B. The redshift of the absorber 
is z = 0.748, consistent with the color of Gl. Therefore the absorber may be associated with the 
lensing galaxy Gl. 
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4. Lens Modeling 

Although the lens systems appear quite complex, we will first attempt to model the lens 
systems with simple models. Specifically, we try models that describe Gl by a Singular Isothermal 
Ellipsoid (SIE), or Singular Isothermal Sphere (SIS) plus external shear. Even these simple models, 
however, have eight parameters (the galaxy position x g and y g , the Einstein ring radius Re, the 
cllipticity e or shear 7, the position angle 6 e or 7 , the source position x s and y s , and the flux of 
the quasar /) that is equal to the number of observational constraints from the UH88 imaging data 
(the image positions, the galaxy position, and fluxes of the images; see Table 2). Thus there are no 
degrees of freedom and in usual cases we will be able to find models that perfectly reproduce the 
observables. To fit the models, we use standard lens modeling techniques as implemented in the 
lensmodel software (Keeton 2001). 

SDSS J1001+5027 is well fitted by both the SIS plus shear and SIE models. The resulting 
fitting parameters are shown in Table 3. The ellipticity e = 0.25 in the SIE model is similar to that 
of the light, though the position angle 8 e = 10.9 (measured East of North) is quite different from 
that of the light (~ —60°) measured from the UH88 image. In general the position angles of the light 
and lens models are aligned (Keeton, Kochanek, & Falco 1998); therefore this result suggests that 
the external field, rather than the galaxy Gl, is responsible for the quadrupole moment of the lens 
potential. The position angles in the models are rather close to the direction to G2. Since in general 
the position angle of external shear gives an idea of the direction to a main perturber, it seems 
that G2 significantly affects the lens models. We also predict the time delay as At ~ 45/i -1 day (A 
leads B), assuming a lens redshift of z = 0.3. 

Fitting SDSS J1206+4332 by either SIS plus shear or SIE failed; the models yielded large 
chi-squares x 2 > 2 with no degrees of freedom, and the models required unnaturally large e or 7. 
We also tried a SIE plus shear model, but the resulting fit was similarly poor. This implies that the 
lens system is too complicated to be described by such simple models. Thus we add G2, which is 
modeled by a SIS, as well as galaxy Gl modeled by SIE, in order to make the model more realistic. 
We derive best-fit models by changing the value of the Einstein radius of G2, Re{G2). We find 
that this "SIE+G2" model fits the data well when R~e(G2) < 1", although the maximum value is 
slightly smaller than that inferred from the Faber-Jackson relation (i?E(G2)/i?E(Gl) = 0.92 from 
the .R-band flux ratio, assuming Gl and G2 are at the same redshift). As a specific example, we 
show the best fit parameters for Re(G2) = 1" in Table 3. The position angle of galaxy Gl in 
the model (9 e = —89.3) is in good agreement with that observed. The time delay is derived to 
be At = 92.6/i _1 day (A leads B), assuming the strong Mg II absorption system at z = 0.748 is 
associated with the lens galaxy. Even if we decrease the value of Re{G2) up to (X'l, the time delay 
is affected only moderately; the time delay is predicted to be 104.4/i~ 1 day when Re{G2) = Qf.'l. 

In summary, lens modeling has revealed that neither system is simple. In particular, the 
secondary galaxies may play an important role in both lenses. 
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5. Lens Galaxy Environments 

Lens galaxies of lensed quasars, particularly in systems with relatively large image separations 
(> 3"), commonly lie in groups or clusters (Keeton et al. 2000; Faure et al. 2004). Such "compound" 
lens systems include Q 0957+561 (Walsh et al. 1979), PG 1115+080 (Weymann et al. 1980), MG 
2016+112 (Lawrence et al. 1984), RX J0911+0551 (Bade et al. 1997), MG 0751+2716 (Tonry 
& Kochanek 1999), SDSS J0903+5028 (Johnston et al. 2003), CLASS B1608+656 (Fassnacht et 
al. 2004), B2108+213 (McKean et al. 2005), and HE 0435-1223 (Morgan et al. 2005). If dense 
environments of the lens galaxies are common, then they could affect strong lens studies in several 
ways (Keeton Sz Zabludoff 2004). On the theoretical side, estimation of environmental effects 
remains controversial: While Keeton et al. (2000) and Holder & Schechter (2003) argued that the 
large fraction of lens systems should lie in dense environments and thus the environmental effects 
are significant, Dalai & Watson (2005) estimated using a halo occupation distribution that the 
typical values for the external convergence and shear are quite small. Although we have already 
seen in §3 and §4 that SDSS J1001+5027 and SDSS J1206+4332 have additional galaxies that likely 
affect the lens potential, it is important to check the larger field for hints of groups or clusters. 

Figure 6 shows the fields around the lens systems obtained at UH88. We find many galaxies 
around SDSS J1001+5027, indicating that there may be a group or cluster along the line of sight. 
For SDSS J1206+4332, we cannot see any noticeable enhancement of the number of galaxies around 
the lens system, though there are several faint galaxies near the lens system. 

To explore the environments further, we derived the number densities of galaxies around each 
lens system using the UH88 images. We perform object identifications using the Source Extractor 
algorithm (SExtractor; Bcrtin & Arnouts 1996). We define galaxies as objects with SExtractor 
parameter CLASS_STAR smaller than 0.6 in the / band image. Note that this star/galaxy separation 
criterion is successful only for objects with / < 22. We derive the number densities of galaxies 
as a function of /-band magnitudes of the galaxies, which are shown in Figure 7. In the Figure, 
we compare the galaxy number densities within 60" of each lens system with those of the back- 
ground, estimated from galaxies more than 60" away in each image after excluding the region near 
bright stars. The total areas of the backgrounds are thus 87.0 arcmin 2 and 79.8 arcmin 2 for SDSS 
J1001+5027 and SDSS J1206+4332, respectively. We find an enhancement of galaxy number den- 
sities for SDSS J1001+5027; the number densities of three magnitude bins show excesses by more 
than la. On the other hand, the galaxy number count around SDSS J1226+4332 is consistent with 
that of the background. 

Figure 8 shows a color-magnitude diagram of galaxies near SDSS J1001+5027 to study the 
origin of the enhancement of galaxy number densities. We find a weak ridge line at R — I ~ 0.75, 
suggesting the existence of a group or a cluster at z ~ 0.2, which is consistent with that of lens 
galaxies Gl and G2 estimated from their colors. We also plot the distribution of galaxies at z > 0.2, 
by making color cuts; V — R > 0.9, R — I > 0.7, and I < 22. We find that the galaxies are clustered 
North and West of the lens system, rather than distributed homogeneously (see also Figure 6). 
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Thus the lens galaxy may lie in a group or a cluster that is located to the north-west of the lens; 
spectroscopic identifications of these galaxies as well as the lens galaxy should be undertaken. 

6. Is the number of 3" lenses consistent with theory? 

Thus far, 13 gravitationally lensed quasars have been discovered using the SDSS. In addition, 
we have recovered several previously known lensed quasar systems. Table 4 summarizes the current 
status of our lens search in the SDSS. Note that limitations of follow-up time have forced us to 
focus on the SDSS spectroscopic sample of quasars 3 . Thus the list does not contain lenses that 
do not have SDSS spectra (e.g., APM 08279+5255; see Pindor et al. 2003). The Table does not 
contain SDSS J1402+6321 (Bolton et al. 2005) either, since the redshift (z = 0.48) is below our 
criterion for the quasar sample (see discussion below). Among gravitational lenses in the Table, 
SDSS J1004+4112 and Q 0957+561 were selected by searching around each quasar for stellar 
objects that have similar colors as the quasar itself (Oguri et al. 2004a; Hennawi 2004), and the 
rest were successfully selected by our standard candidate selection algorithm (see §2) though the 
first identifications of some of the candidates were made by different algorithms. An exception is 
SDSS J0903+5028 (Johnston et al. 2003) which was targeted as a luminous red galaxy (Eisenstein 
et al. 2001) because of the bright lens galaxy; however, we include this system in the statistical 
analysis below, since the quasar would have been targeted as a high-redshift quasar if it had not 
been obscured by the foreground lens galaxy (Richards et al. 2002). The Table shows that we 
have discovered a relatively large number of ~ 3" lenses. We note that our current lens sample 
shown in Table 4 is quite incomplete 4 , and that future follow-up observations would increase the 
number of lenses even in the current quasar sample. Therefore we now calculate the expected 
lensing probability in the SDSS quasar sample and compare it with the number of ~ 3" lenses in 
the sample. In particular, we neglect the contribution of lens galaxy environments, as has normally 
been done, to see whether the assumption is still valid or not. In computing the lensing probability, 
we assume a spatially flat universe (Qm + = 1). 

We compute the lensing probability distribution along the lines described by Turner, Ostriker, 
& Gott (1984). The lens galaxies are modeled by SIS p(r) = a 2 /(2irGr 2 ), where a is the velocity 
dispersion of the lens galaxy. We adopt the velocity function of early-type galaxies determined 
from ~ 30, 000 SDSS galaxies at 0.01 < z < 0.3 (Sheth et al. 2003; Bernardi et al. 2005; Mitchell et 
al. 2005), and neglect the redshift evolution (i.e., we assume that the velocity function is constant 
in comoving units). The use of the velocity function of early- type galaxies is sufficient for our 



The photometric quasar sample (Richards et al. 2004) has an order of magnitude more quasars than the spec- 
troscopic sample, and thus is expected to contain many more lenses. 

4 We note that our follow-up observations are biased against discovering lensed quasars with separations < 2" 
because of poor seeing in many follow-up observations. This mostly explains the large discrepancy between theory 
and observation at the image separation range. 
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calculation because at 9 > 3" lensing by early- type galaxies is dominant (Turner et al. 1984). 
We need also to incorporate the selection function of the SDSS lens search; which we do in a 
preliminary way by making the following assumptions. First, we use the magnification factor of 
the brighter image, ^bright = (#e/#s) + 1) where 9e is the Einstein radius and 9s is the position of 
the source relative to the lens galaxy, to compute the magnification bias, because at 9 > 3" two 
lensed components are well separated in the SDSS data. Although we may have an alternative 
choice of using the magnification factor of the fainter image, we adopt /^bright to calculate the 
upper limit of the lensing probability; thus the actual lensing rate at 9 > 3" might be smaller 
than our calculation 5 . Next, for the limiting bright-to- faint flux ratio we assume / max = 10; this is 
justified because all lenses in Table 4 satisfy this condition. Lensing probabilities calculated with 
this selection function represent an upper limit at relatively large image separation 9 ~ 3". 

To calculate the expected number of lensed quasars in the current SDSS quasar sample, we 
need the luminosity function of quasars as well as the redshift and magnitude distributions of the 
quasar sample. We adopt a sample of ~ 47, 000 quasars at 0.6 < z < 4.0. The sample is constructed 
in the same way as in Oguri et al. (2004a). We did not search for lensed quasars at z < 0.6 as 
low-redshift quasars are intrinsically extended, and we do not use high-redshift quasar {z > 4.0) 
sample because it contains a significant fraction of objects whose spectra were misidentified by the 
SDSS spectroscopic pipeline (c.f., discussion in Schneider et al. 2003). The redshift distribution is 
similar to that shown in Oguri et al. (2004a). For the luminosity function of quasars, we use what 
is called LF1 in Oguri et al. (2004a), which has a faint end slope of 1.64 and a bright end slope of 
3.43 (z < 3) or 2.58 {z > 3). 

The result is shown in Figure 9. As seen, the number of ~ 3" lenses already exceeds the 
theoretical expectations, when we adopt the standard value of the cosmological constant, Q\ ~ 0.7. 
The situation is similar if we increase f^A to 0.8. In terms of the constraint on the cosmological 
constant, the excess in the bin centered at 9 = 3" is described as JIa > 0.90 (68% confidence 
limit) assuming a Poisson distribution, which is highly inconsistent with recent measurements of 
the cosmological constant Q\ ~ 0.70 with errors of less than 10% (e.g., Tegmark et al. 2004). We 
emphasis that our current lens sample is quite incomplete; indeed, currently we have several ~ 3" 
quasar pairs with the same redshifts which could also be gravitational lenses (Hennawi 2004), which 
need deep imaging to find the putative lens galaxies. Thus observations of these quasar pairs as 
well as lens candidates should be conducted, which might make the discrepancy even larger. Such 
an excess may indicate that our simple treatment of lensing statistics is inaccurate and that we 
need to include other effects. 



5 Note that at smaller separation (8 ~ 1") our computation is expected to underestimate the lensing probability 
because at such small separations lensed components are not decomposed and thus the magnification factor of the 
two images pitotai = 2#e/#s is more appropriate. 
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7. Summary 

We report the discovery of two gravitationally lensed quasars, SDSS J1001+5027 and SDSS 
J1206+4332. The systems were identified as new lens candidates in the SDSS, and confirmed as 
lenses by spectroscopic and imaging observations at the ARC 3.5m and the UH88 telescope. SDSS 
J1001+5027 is a lensed quasar at z = 1.838, and consists of two lensed images separated by 2'.'86. 
SDSS J1206+4332 is a lensed quasar at z = 1.789, and consists of two lensed images separated by 
2'.'90. In each system we have identified the galaxy responsible for the lensing. 

We have found that the lens systems are complicated. The imaging data clearly show other 
galaxies that are very close to the main lens galaxies. The lens modeling has shown that these 
galaxies affect the lens potentials significantly. We have examined the wide field images, which 
show an enhancement of the galaxy number density within 60" of SDSS J1001+5027. Spectroscopic 
follow-up observations are needed to see if the group/cluster (z ~ 0.2) is actually associated with 
the lens galaxy. 

Although the SDSS lens survey is ongoing, we have made a preliminary comparison of the- 
oretical lensing probability distributions with the observed distribution. We have found that the 
number of lenses with 6 ~ 3" already exceeds the theoretical expectations. We still have many 
lens candidates with 9 > 3" that remain to be observed (Hennawi 2004); if some of them turn 
out to be true gravitational lenses, the conflict will become even stronger. This excess may be 
caused by external convergence and shear fields which we have not taken into account in our cal- 
culation. Basically, external shear broadens the distribution of image separations for a given mass 
of a lens object (e.g., Huterer, Keeton, & Ma 2005). This broadening is enough to enhance the 
lensing probability at 9 > 3" because the lensing probability at the image separation region is a 
strong function of image separations. More significant enhancement may be achieved by external 
convergence, since it increases both image separation and lensing probability. Indeed, among five 
intermediate-separation lenses in Table 4, two lens systems (RX J0911+0551 and Q 0957+561) lie 
in clusters, and two other lens systems (SDSS J0903+5028 and SDSS J1001+5027) also appear to 
lie in dense environments. In addition, the two lens systems reported in this paper are complex in 
the sense that galaxies very close to the main lens galaxies affect the lens potentials. Other possible 
systematic effects include triaxiality of lenses (Oguri & Keeton 2004) and massive substructures 
(Cohn & Kochanek 2004). In either case, lens statistics at this image separation should be done 
with caution; simple models that consider only isolated single lens objects can be misleading. 
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Table 1. SDSS Astrometry and Photometry 



Object 


R.A. (J2000.0) 


Dec. (J2000.0) 


u g 


r 


i 




SDSS J1001+5027 


A 


10 01 28.61 


+50 27 56.9 


17.71 + 0.02 17.60 + 0.02 


17.55 + 0.03 


17.36 ± 0.03 


17.33 ± 0.04 


B 


10 01 28.35 


+50 27 58.5 


18.60 + 0.08 18.35 + 0.08 


18.11 + 0.06 


17.71 ± 0.06 


17.62 ± 0.05 


SDSS J1206+4332 


A 


12 06 29.65 


43 32 17.6 


18.66 + 0.03 18.80 + 0.05 


18.75 + 0.05 


18.53 ± 0.04 


18.48 ± 0.04 


B 


12 06 29.65 


43 32 20.6 


19.69 + 0.08 19.43 + 0.12 


19.33 ± 0.08 


19.22 + 0.08 


18.83 + 0.04 



Note. - PSF magnitudes returned by the SDSS photometric pipeline are shown. Reddening corrections are not 
applied. 
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Table 2. UH88 Astrometry and Photometry 



Object 


x[arcsec] a 


y[arcsec] a 








SDSS J1001+5027 


A 


0.000 ± 0.005 


0.000 ± 0.005 


18.12 ±0.01 


17.69 ±0.01 


17.32 ±0.01 


B 


2.418 ±0.009 


1.526 ± 0.005 


18.53 ±0.01 


17.98 ±0.01 


17.48 ±0.01 


Gl 


1.779 ± 0.049 


0.857 ±0.123 




20.51 ±0.03 


19.63 ±0.03 


G2 


1.795 ± 0.088 


-0.700 ± 0.053 




20.91 ±0.04 


20.15 ±0.03 


G1+G2 C 






22.09 ±0.08 






SDSS J1206+4332 


A 


0.000 ±0.011 


0.000 ± 0.010 


18.63 ±0.02 


18.54 ±0.01 


18.05 ±0.02 


B 


-0.098 ± 0.006 


2.894 ± 0.009 


19.12 ±0.02 


18.95 ±0.02 


18.38 ±0.02 


Gl 


-0.664 ±0.137 


1.748 ± 0.028 




21.34 ±0.07 


19.51 ±0.03 


G2 


1.320 ±0.147 


5.999 ±0.148 


22.43 ±0.13 


21.16 ±0.08 


19.63 ± 0.04 


G3 


-2.052 ± 0.200 


2.397 ±0.152 


21.78 ±0.09 


22.05 ±0.10 


22.57 ±0.11 



a Positions relative to the brighter quasar components are presented. The positive di- 
rections of x and y are defined by West and North, respectively. Error bars do not include 
the error of the image scale. 

b Error bars do not include the zero-point error of ~ 0.1 mag. 

c In the F-band image it was difficult to decompose Gl and G2 because of the faintness 
if the image, thus we show the total magnitude of these two galaxies. 
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Table 3. Lens modeling 




Model 


i?E[arcsec] e or 7 6 e or # 7 [deg] a 


A*[/r 1 day] b 


SDSS J1001+5027 


SIS+shear 
SIE 


1.35 0.09 11.9 
1.38 0.25 10.9 


42.4 
47.0 


SDSS J1206+4332 


SIE+G2 C 


1.39 0.31 -89.3 


92.6 



Note. — See §4 for details of the lens models adopted here. 

a Each position angle is measured East of North. 

b For the redshifts of lens galaxies, we assumed a lens redshift 
of z = 0.3 for SDSS J1001+5027, and the redshift of strong Mg II 
absorption, z = 0.748, for SDSS J1206+4332. We also assumed 
Om = 0.3 and = 0.7, but the values are quite insensitive to 
these cosmological parameters. 

c The Einstein radius of G2, i? E (G2), is fixed to l". 
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Table 4. Gravitationally Lensed Quasars (Re-)Discovered in the SDSS: Current Status 



Name 


f\ 'A 

u max 


Z s 


Zl 


Comments 


Rcf. 


Keg ui^tz — 1UU 


9" 9 9 
z. zz 


9 79 
Z. 1 z 


n aq 


picvioiiDiy Known 16110 


1 
1 


C!r>C!C! TD9zLfi 089^ 


1"D/L 

1 . U^i 


i.oy 


n 79? 

U. ( Z : 




9 
z 


DiJoj JUyuo-rOUZo 


9" 8 3 
Z. OO 


3 ^8 
o.Oo 


u.oy 


IdlgeXcCL do ex gdicxXy 


9 



l\jY JUyil-|-U001 


O. zo 


9 an 

Z.oU 


n 77 


previously Known lens 


zl 
4 


odd uyuy-t-ozo 


1"1 1 

1.11 


1 38 

1.00 


n 83 


pieviousiy Known lens 


K 



huhb J09z4+oziy 


1. 78 


1.5z 









Q 0957+561 


6'/ 17 


1.41 


0.36 


previously known lens 


7 


SDSS J1001+5027 


2?86 


1.84 






8 


SDSS J1004+4112 


14'.'62 


1.73 


0.68 


large-separation lens 


9 


SDSS J1021+4913 


l'/05 


1.72 






10 


PG 1115+080 


2'/43 


1.72 


0.31 


previously known lens 


11 


SDSS J1138+0314 


l'/46 


2.44 






12 


SDSS J1155+6346 


l'/83 


2.89 


0.18? 




13 


SDSS J1206+4332 


2 / /90 


1.79 


0.75? 




8 


SDSS J1226-0006 


l'/24 


1.13 


0.52? 




14 


SDSS J1335+0118 


l'/56 


1.57 






15 


SDSS J1650+4251 


l'/18 


1.55 






16 



Note. — The SDSS strongly lensed quasar survey is still ongoing, and we 
have many other candidates requiring follow-up observations. Therefore, 
this list is far from complete. This list includes only lenses that have SDSS 
spectra. 

a The maximum separation between multiple images. 

References. — (1) Surdej et al. 1987; (2) Inada et al. 2005b; (3) Johnston 
et al. 2003; (4) Bade et al. 1997; (5) Oscoz et al. 1997; (6) Inada et al. 2003a; 
(7) Walsh et al. 1979; (8) this paper; (9) Inada et al. 2003b; (10) Pindor 
et al. 2005; (11) Weymann et al. 1980; (12) Buries et al. 2005; (13) Pindor 
et al. 2004; (14) Inada et al. 2005a; (15) Oguri et al. 2004b; (16) Morgan, 
Snyder, & Reens 2003 
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Fig. 1.— The SDSS i-band images of SDSS J1001+5027 (left) and SDSS J1206+4332 (right). The 
image scale is CK'396 pixel" 1 . 
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Fig. 2. — ARC 3.5-m spectra of SDSS J1001+5027 components A and B. The exposure time was 900 
seconds. Both components have the same redshift, z = 1.838. The strong absorption at ~ 7600 A 
is atmospheric. The feature at ~ 6200 A is a bad column. The spectrum is smoothed by a 3 pixel 
boxcar. The ratio of the spectra as a function of wavelength is shown in the bottom panel. B is 
evidently reddened. 
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UH88/I 
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•* 


A 


1 

L . 












* 

" * 


G2 
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Fig. 3. — UH88 and ARC 3.5-m images of SDSS J1001+5027. Lower panels are same as upper 
panels, but after subtracting two point source quasar components. The image scales are 0'.'232 
pixel" 1 (UH88) and 0"282 pixel" 1 (ARC), and the seeings were ~ 0"8 (UH88) and ~ l'.'O (ARC). 
North is up and East is left in all panels. We find two galaxies in the residual images shown in the 
lower panels, which are labeled Gl and G2. 
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Fig. 4. — ARC 3.5-m spectra of components A and B of SDSS J1206+4332. The exposure time was 
1800 seconds. The redshifts of both components are the same, z = 1.789. The strong absorption 
at ~ 7600 A is atmospheric. The spectrum is smoothed by a 3 pixel boxcar. There is strong Mg II 
absorption at ~ 4900 A (corresponds to an absorber at z = 0.748) in the spectrum of component 
B. The ratio of the spectra as a function of wavelength is shown in the bottom panel. 
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Fig. 5. — UH88 images of SDSS J1206+4332. Lower panels show images after subtracting the two 
quasar components. North is up and East is left in all panels. We find three galaxies in the residual 
images shown in the lower panels, which are labeled Gl, G2, and G3. 
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Fig. 6.— The 2' x 2' /-band fields around the lens system SDSS J1001+5027 (left) and 
SDSS J1206+4332 (right) taken at UH88. The exposure times were both 360 seconds. North 
is up and East is left in both panels. 
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Fig. 7. — /-band galaxy number counts of fields within 60" of the lens systems (filled triangles with 
solid errorbars) as well as background counts (crosses with dotted errorbars). The errorbars include 
Poisson errors only. The vertical dashed lines indicate the approximate limit of the star/galaxy 
separation. 
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Fig. 8. — Left: Color-magnitude diagrams of galaxies in the field within 60" of SDSS J1001+5027 
(filled triangles) and the background (crosses). Right: The distribution of galaxies at z > 0.2 (filled 
triangles) around the lens system SDSS J1001+5027 (open square). North is up and East is left. 
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Fig. 9. — Tentative comparison of image separation distributions of lensed quasars in the current 
SDSS quasar sample. Lines show the theoretical expectations (see text for details of computation) 
using an SIS model, for different values of the cosmological constant, Qa = 0.6, 0.7, and 0.8 from 
bottom to top, respectively. Shaded region indicates the distribution of all lensed quasars discovered 
and recovered in the SDSS (see Table 4). The large-separation lensed quasar SDSS J1004+4112 
(6 = 14'.'6) is not included in this plot. 



